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Abstract 

The quantum mechanical version of a classical model for studying the orientational 
degrees of freedom corresponding to a nematic liquid composed of biaxial molecules 
is presented. The effective degrees of freedom are described by operators carrying an 
SU(3) representation, which allows the explicit calculation of the partition function 
in the mean field approximation. The algebraic consistency conditions are solved nu- 
merically and the equilibrium phases of the system are determined. In particular, the 
entropy, the specific heat and the order parameters are presented for different choices 
of the constituent biaxial molecules. Our results reproduce the classical calculation in 
the limit of high temperatures and high quantum numbers. 
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1 Introduction 



Nematic liquids constitute a particular case of the phases denominated with the generic 
name of liquid crystals 0, |2|- This is a state of matter which is intermediate between the 
liquid and the solid phases. It lacks the positional symmetry characterizing the lattice of a 
crystal, which makes it similar to a liquid, but it retains some orientational symmetry that 
is not present in a normal (isotropic) liquid. The molecules forming a liquid crystal can be 
roughly caracterized as rod-like or disc-like, and the corresponding orientational symmetry is 
manifested by different types of local aligments of such molecules with respect to a preferred 
direction characterized by a unit vector called the director. In the case of a nematic liquid, 
the director is fixed in space and the molecules remain, on the average, parallel to it, though 
they can difuse as in a normal liquid. The existence of a preferred direction makes this phase 
to exhibit birrefringent properties. On the other hand, the energies involved in the transition 
from the liquid crystal phase to the isotropic liquid phase are normally small compared with 
the corresponding ones in the solid-liquid transition. This makes it possible to drastically 
change the response of liquid crystals by external agents, using relatively small amounts 
of energy. It is this property which determines the numerous and important technological 
applications of liquid crystals 0]. 

Besides the important applications of liquid crystals in everyday's life, together with 
the huge amount of theoretical work that has been consequently developped, it has been 
suggested recently that these phases might be relevant also in other situations, like for 
example, the description of nuclear matter under very stringent conditions of densities and 
temperatures [|J, particularly in the case of neutron stars. These quantum objects are cold 
stars composed primarily of neutrons and supported against collapse by neutron degeneracy 
pressure. In this situation, it is predicted that the inner crust of a neutron star (roughly 
located at 9.2 km from the center of a typical neutron star with a radious of 10 km and a mass 
of 1.5 solar masses) is separated from the core, composed mainly of uniform nuclear matter, 
by a thin spherical shell about 0.1 km wide, containing nuclear matter in the form of non- 
spherical nuclei (rods, plates, tubes) plus dripped neutrons plus electrons. This liquid crystal 
behavior, as opposed to a crystalline solid, will have important consequences for a number 
of aspects of neutron star properties ||. Nevertheless, there is a main difference between 
this case and the situation in the laboratory: in the neutron star case, the basic entities 
from which nuclei are constructed are nucleons, which are spherical; while in laboratory 
liquid crystals, the basic ingredientes are non-spherical biaxial molecules in general ||. In 
this respect, it might be interesting to consider liquid crystals made up from constituentes 
having permanent magnetic moments, as nucleons, generating a liquid ferromagnet. In the 
laboratory, such structures, though not forbidden, should be very rare because the exchange 
couplings beteween neighboring magnetic moments are much weaker than the conventional 
interactions between molecules. 

Motivated by some of the above considerations, we consider the quantum mechanical 
description of the orientational degrees of freedom of a nematic liquid consisting of biaxial 
molecules, in the mean field approximation. This system has been previously discussed in 
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Ref. 0, where an exact calculation of the classical partition function, together with an exten- 
sive analysis of the corresponding phases, transition temperatures, themodynamic quantities 
and order parameters, was made. We will refer to this work as the classical model. In this 
reference, the effective orientational degrees of freedom of the non-spherically symmetric 
molecules were described by an SU (3) invariant Hamiltonian constructed from the angular 
momentum operators Lj, i = 1,2,3 together with the traceless quadrupole tensor operator 
Qap, a,P = 1,2,3, which constitute a representation for the generators of 577(3). The 
interpretation of the model was performed in a coordinate system where the molecule is at 
rest, i.e. Lj = 0. In this system, the quadrupole operator Q a p relates to the inertia tensor 
of the molecule in such a way that the 577(3) Casimir operators, related to the eigenvalues 
°f Qa/3, provide a description of the shape of the molecule. 

The classical partition function in the mean field approximation was exactly calculated 
in Ref.JJ by first rewritting it as an integral over the space parameters of the 577(3) group, 
according to a theorem due to Macfarlane et al. 0. This integral has been previously 
calculated in the bosonic case by Harish-Chandra, and Itzykson and Zuber ||. The exact 
evaluation of the partition function, together with the fact that the consistency conditions 
consisted in two coupled algebraic equations, constituted an advantage over previous sta- 
tistical theories proposed |J in order to explain the isotropic — > uniaxial — > biaxial phase 
transitions discovered experimentally in nematic liquids |10[| . The works included in || are 
generalizations of the Maier-Saupe theory, that was for some time the fullest existing statis- 
tical mechanical theory of the nematic state, which, nevertheless, did not include the biaxial 
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Following an approach similar to some previous works, we only consider the orientational 
interaction energy among the molecules of the liquid, in such a way that the exact nature 
of the intermolecular forces need not be specified . Whether the liquid crystal state arises 
because of the hard core repulsive forces or the long range dispersive forces is a matter 
of indifference in the following. This approximation is expected to produce qualitatively 
reasonable results in the nematic case ||. In this way, we are really calculating the excess 
thermodynamical properties of the ordered system, relative to those of the disordered one. 

In this work we present a quantum description of the problem discussed in [^|, which 
amounts to the generalization of the usual Heisenberg model describing the interaction of 
angular momentum operators (SU(2) case ) to the situation where the interacting operators 
carry an 5?7(3) representation. Similarly to the classical case, we are able to exactly calculate 
the partition function in the mean field approximation. This is achieved by using the Weyl 
formula for the characters of 577(3) [Tj|. We expect to recover the results of Ref. in the 



high temperature, high quantum numbers limit [jl3| , |14 . 

The paper is organized as follows. In section 2 we introduce the model and apply the 
Weyl formula to calculate the partition function. Section 3 describes the characterization of 
the constituents molecules of the liquid in terms of the irreducible representations of SU(3). 
The equilibrium phase is parametrized in the same way as in Ref. as briefly recalled in 
section 4. The connection with the classical case is presented in section 5, where a function 
relating the corresponding partitions functions is explicitly calculated. Section 6 contains 
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the selfconsistent equations that allow for a numerical determination of the equilibrium 
parameters, which are obtained by determining the absolute mimima of the free energy of the 
system. Besides, some symmetry properties of the quantum partition function are exhibited 
there. Section 7 describes the calculation of the equilibrium parameters of the different phases 
of the system, emphasizing those limiting cases that can be treated algebraically. We have 
also calculated numerically the specific heat and the entropy as functions of the adimensional 
temperature, for some representative cases. The order parameters in the quantum case are 
also defined in this section. In many of the figures presented, we have shown both the 
classical and quantum predictions in order to assess the main differences between the two 
approaches. Some conclusions and comments are presented in section 9. Finally, there is 
an Appendix summarizing some basic properties of SU(3). Here we also give an alternative 
way of rewritting the partition function, in terms of polynomials, which have been useful in 
obtaining some limiting cases considered along this work. 



2 The model 

Following the work of Ref . || , let us consider the effective Hamiltonian 

h = -^Y, sts:, (i) 

which describes the orientational energy of the system. Here, the N independent operators 
S™, m — 1, 2, 3, . . . , N, a = 1,2,3, ... ,8 describe the degrees of freedom of each molecule of 
the liquid and carry an SU(3) representation, i.e. 

[S™, S£] = 8 mn if abc S?, (2) 

where f abc are the SU(3) structure constants. As it is well known, the fundamental repre- 
sentation of this group is written in terms of the Gell-Mann 3x3 matrices A a given in the 
Appendix, where other useful properties are also recalled. 

In the mean field (molecular field) aproximation we replace ([!]) by 

H = -JB a Y / S n a , (3) 

n 

where the eight-dimensional vector B a = ■h < J2 m >, which represents the average field 
seen by each molecule, has to be determined consistently. 
In this approximation, the partition function Z is 

Z = Tre-? H = Tre? JB « ^ s " = Z» ' (4) 

where 

Z = Tre pjBaSa . (5) 
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As usual -p = kT, where k is the Boltzman constant and T is the temperature of the system. 



The consistency equations for the mean field approximation reduce to 

1 d 



B„ 



(3JdB a 



-InZn 



(6) 



The basic problem now is to calculate such partition function for an arbitrary (p, q) 
irreducible representation of SU(3). The first step is to select and adequate coordinate 
system to perform the calculation. Since the trace is invariant under a unitary transformation 
of the operator inside, we can diagonalize the hermitian operator j3JB a S a in such a way that 
it can be written as (3J{B 3 S 3 + B 8 Ss), where we assume that the operators S3 and Sg are 
chosen to be simultaneously diagonal. This is equivalent to choose a coordinate system where 
the eight-dimensional vector B a has only non-zero components in the 3 and 8 axis. The fact 
that we need two non-zero components, instead of one as in the 577(2) case, is related to 
the property that in the adjoint representation of SU(3) there are two invariants: B a B a and 
d a b c B a B b B c . Using the explicit form of d a f, c , it can be shown that it is possible to select the 
B 3 and B 8 components as the only non-zero ones. In this way we are left with 



Z = T re PJ{B 3 S 3 +B 8 S 8 ) _ 

From now on we will assume that B a = (0, 0, B 3 , 0, . . . , B 8 ). 

In the fundamental representation (1,0), the partition function is given by 



(7) 
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JJE-jU 



(8) 



where /2j are the weights of the fundamental representation. 

The partition function @ is the character of the operator e~^ H . In this way, to calculate 
Z Q for an arbitrary representation (p, q) of SU(3), we can use the Weyl formula for the 
characters of SU(n), given in terms of the eigenvalues of the operator in the fundamental 



representation. For SU(3), the Weyl formula reduces to [12 



Yj{»?i ,m} 



(A x -A 2 ) (Ax - A 3 ) (A 2 -A 3 





Af +1 
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Ar 2 
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A^ 2 


Af +1 


1 



(9) 



Here rji = p + q is the number of boxes in the first row of the corresponding Young 
tableau, while — Q is the number of boxes in the second row of the tableau. In other 
words, p is the number of columns with one box (p > 0) and q is the number of columns 
with two boxes (q > 0) in the tableau. The representation (0, 0) is the trivial one. The 
eigenvalues Aj, i — 1,2,3 are given by 

A . = e PJS.Hi ( 10 ) 
and satisfy the condition A!A 2 A 3 = 1. The weights /2j are given in the Appendix. 
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A more convenient notation is 

P=\PJB 3 , a = ^=pjB 8 , (11) 

which produces the following expressions for the eigenvalues of the Hamiltonian in the fun- 
damental representation 

A 1 = ef+°, A 2 = e~ p+a , A^ 1 = A,A 2 = e 2 ° . (12) 

The partition function is given by 

4 P ' 9) = — r~ 9 ^ (-e^^ sinh (p + l)p 

sinn 2p — 2 smh p cosh 6a v 

+e {p - q)a smh{p + q + 2)p-e- {q+2p+3)a sinh {q + l)p) . 

(13) 

It is convenient to observe that the above expression for zjf' q ^ can be rewritten as 

^ = iw> <"> 

where 

N (P,Q) = e (p-q)c sinh ( p + q + 2 )p - e (P+29+3)- sinh ^ + ^ p 

- e -(9+ 2 P+3)- S inh(g+ (15) 

3 Specification of the constituent molecules 

The orientational degrees of freedom of the molecules are described by the eight operators 
S a , which can be reinterpreted as the angular momentum Li, % — 1, 2, 3 and the quadrupole 
moment Qij operators, according to ||, 

L\ = 25V, L 2 = —25*5 L3 = 25*2, 
Q12 — Si, Q13 = S4, Q23 = Sq, 



Qll — S3 + ~^Ss, Q22 — —S3 + ~^S8, Q33 — ~~~^^ 8 ' 



(16) 



The molecules forming the liquid will be characterized by the representation (p, q) of £77(3). 
To this end we consider the two Casimir operators of the group 

C2 = S a S a = J2I, C3 = d a b c S a SbS c = I3I (17) 
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which take the following values in each irreducible representation (p,q), |T5| 

p 2 + q 2 + pq + 3 (p + q) 



;is) 



h = ^(p - g)(p + 2g + 3)(2p + g + 3). (19) 

Motivated by the classical description and in order to compare with it, we introduce 
a parametrization of the above invariants in terms of fa, % = 1,2,3, with the restriction 
/xi + /i2 + to = 0. In the classical case one defines 1$ = ^tototoi which can be reproduced 
in our case by choosing 

Hi = ~(p + 2q + 3), H2 = ~(j>-q), to = \{^P + q + 3), (20) 



according to the relation (19) . Following the analogy, we next calculate what used to be the 
quadratic Casimir operator. We obtain 

- (toto + toto + toto) = h + 1, (21) 

where I2 is given by Eq. ([18D . 

The choice fl2~0D implies the ordering fa < /i 2 < A*3> f° r arbitrary values of p and q. 
The eigenvalues (|20"D satisfy the cubic equation 

fj*-(l 2 + l)ti-ll 3 = 0, (22) 

and can be parametrized as 

fit = ix cos + , to = fx cos (if) - ^pj , to — ^ cos ^ (23) 

which is consistent with the ordering introduced in (p0|) . To this end we must always have 

\x sini/j > 0, [i cosif) > 0, tamp < v3, (24) 

where 

/ (/ 2 + l) , y/3/3 , ,, 

V = ±2 \ O ' COS3 ^ = ± 3"' ( 25 ) 

V 3 (h + 1) 5 

in complete analogy with the classical case. The above conditions lead to the following 
correlations between if) and \i 

/x> <ifj < -, Li<0<-+7r<ip<7r + -, (26) 

for arbitrary p and g. In this way, for p and q fixed, the change fx — > — /x has to be followed 
by -0 — > ^ + 7r. Thus we pass from the upper plus sign to the lower minus sign in the two 
relations contained in fl25|). 
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The physical interpretation associated to the (p, q) characterization of the individual 
molecules is made through the classical connection 



1 1 

a* = -(at 2 + a 2 2 + a 3 2 ) + (27) 



which relates the /x-values ( |2"0"D with the shape corresponding to a molecule described by a 
mass distribution defined by six equal masses at distances ±ai, ±a2, ±03 from the center [[J. 
Here we are assuming that the a^'s are written in dimensionless units. The correspondence 
is 

A*i = A*2 = A*3 (spherically symmetric molecule), 
j"i = A*2 < 0, /i3 > 0, (rod molecule), 
/ii < 0, /i 2 = /i 3 > 0, (disc molecule), 

-(//1 + /i 3 ) > // 2 , (rod— like molecule), 

+ ^3) < A*2> (disc— like molecule). (28) 



The expressions (20) together with the classification (28), imply that in the quantum 
case we do not have at our disposal either spherically symmetric or uniaxial molecules to 
begin with. The description used allows only biaxial molecules which are rod-like (p > q) or 
disc- like (p < q). 

Using the parametrization (p3| ) we obtain the following expressions for the functions of 
p, q appearing in the partition function ([13]) 



1 2n 1 27T 

-g(P - Q) = V> C0S W - y ), VP + g + 2) = /i sm(^ - — ), 

1 27T 1 27T. 

--(p + 2g + 3) = fi cos{ip + —), -j={p + 1) = fi sm(ip + —), 

+ ^(2p + q + 3) = fi cosip, -^=(q+ 1) = fj, sin-0. (29) 
3 a/3 



4 Parametrization of the equilibrium phases 

The equilibrium state of the system is characterized by the vector B a in the frame where its 
only components are B 3 and B s . Again, an SU(3) invariant description of such phase will 
be given by the corresponding vector invariants, in the adjoint representation 

h = B a B a = B 3 2 + B$ 2 , I 3 = d a b c B a B b B c = -j=(—B$ 3 + 3B 8 B 3 2 ). (30) 

These invariants take continuous values as oppossed to those describing the individual 
molecules, which have only discrete values. 



7 



From Eqs (0), we can verify that Bg satisfies the cubic equation 

2, 



B 8 3 - I 2 B 8 - |/ 3 = (31) 



Ml • 


■B = 


Ai = 




■ B = 


A 2 = 


#3 ' 


B = 


A 3 = 



and consequently can be parametrized in any of the forms (|2"3"D. A convenient choice is made 
by considering that the relevant variable is the matrix B a X a in the fundamental representa- 
tion, which eigenvalues are 

1 1 / 2tt 

1 1 / 2tt 

"2 S3 + 27! s:=Acos r~T 

— L^g ;= \ cos 0. (32) 
V3 

In this way we identify 

B 8 = —V3 A cos 0, B 3 = —V3 A sin0. (33) 

The equilibrium phase is described by the variables A, <p, which can be translated in a ge- 
ometrical language regarding the average shape of the equilibrium states via the relations 
(|28|) and (|32|), in terms of the eigenvalues Aj. 



5 Connection with the classical case 

The parametrization fl29|) allows us to rewrite the quantum partition function in a form that 
closely resembles the classical case. In fact, the factor N^ p ' q ^ of Zo^' g \ given in Eq. fljjD , can 
be rewritten as 

N (r«) = _ V.) - ±A(X, <P + ^)= l -D(X, 0, V), (34) 

where X = ^fiy and D(X,<J>,ip) is the same expression as in the classical case, since we 
have 

A(X, a) = e Xcosa + e Xcos ( a -^ + e Xcos ( Q +¥). (35) 
Here we have introduced the following dimensionless variables 

Jj* y = — x > p = ~2 sm( f ) ^ a = ~273 co ^ 

The temperatute T c used in Ref. || is related to our definition by 

T C =^T, (37) 
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where /i is given by Eq. (|25|) . In the sequel we will refer to T c as the classical temperature. 
The denominator of the partition function can also be written as 

jV(°>°) = -4 sinh(-sin0) sinh(-sin((/) - — )) sinh(-sin(0 + — )), (38) 
2 2 3 2 3 

which allows a more direct comparison with the classical case. 

In fact, in the limit T — > oo (y — > 0, X — > 0) the above expression reduces to 

1 2n 2n 1 

#(0,0) = — y 3 sin( ^ gin ^ j gin ^ + _) = _ y 3 sin3 0_ ( 39 ) 

2 3 3 8 

Since in the case p = q = we have \i = y = X, I 3 = I 2 = 0, sin3ip = 1, we can write 
jy(o,o) _ i sm 30 sinS^. Then, in this limit we have 

A 05 sm30 sinoif) 

In order to compare the full quantum partition function with its classical counter part, we 
subsequently make a precise connection between the two approaches, by rewriting the vari- 
able x c = |A c yu c used in ||, in our terms. Here \i c denotes the "radius" in the parametrization 
of the constituent phase eigenvalues, and fi c = fi q , according to Eqs.pOj). The subindices 
c,q refer to the classical || and quantum cases, respectively. The value A c is the "radius" 
of the parametrization of the equilibrium phase, given by the eigenvalues Aj in Ref 0. The 
variable S a of Ref. || is what we call B a and the matrix whose eigenvalues are the Aj's 
is Mi = £yB a (2\ a ), where we have used that (A a ) c = 2(\ a ) q . We have labeled the "ra- 
dius" associated to the parametrization of the eigenvalues of B a X a by X q = J^^j in such a 

way that A c = (3JX q . Then we obtain x c = |/3JA g /x 3 . On the other hand, we have that 

X = ^fAq y = \f3JX q n q = x c . In this way, fllPf ) is exactly the partition function given in 
Ref §. 

Thus, we can readily see that the numerator of our quantum partition function (|TJ) 
coincides with the classical case, except for the fact that the angle tp is rectricted to discrete 
values according to the relations (|18|), Ql9|) , (p5j), (|20|) and (p3|) . The denominator is quite 
different from the classical case and shows no dependence upon the angle ip, i.e. upon the 
input parameters p, q. 

In general, we can relate the quantum [Z^ 9 ^ ) and classical {{Zq) c = Q) partition 
functions by 

Z = ZQ, (41) 

with 

z = D M fsm(f ~ 4>) fsin(f + 0) 

sinh(|sin0) sinh(|sin(^ L — 0)) sinh(|sin(y : + (/))) ' 

where D^ p ' q ^ is de dimension of the (p,q) representation of 57/(3). Let us observe that the 
quantum calculation depends upon the parameter /x, as oppossed to the classical case where 
this variable appears only through the combination x c . 
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6 Calculation of the equilibrium parameters 



In the coordinates chosen so far, the consistency equations in the mean field approximation 
are given by 

B 3 = —JL\ nZo ^\ B 8 = —JL\ nZo ^\ (43) 
3 pJdB 3 ' 8 f]JdBs K J 

where the partition function Z^^ is written in Eq . (|i~3|) . The original expression 



B a = ^Tr (S a e^ B ^ +B ^) , a = 3, 8, 
z v 7 



(44) 



for the consistency conditions, allows us to see that we will always have the solution B 3 = 
= B s , (A = y = 0) corresponding to the isotropic phase. In this case, the partition 
function is Zq = D^ p ' q \ independent of the variable (f), which we choose as = to describe 
the isotropic phase. Nevertheless, this phase will not always be the stable one. 

Changing to the variables A, (j) introduced in Eq . (j33|) , the equilibrium equations reduce 

to 



A 



1 d 



lnZ (p ' 9) , lnZ (p ' 9) = 



which can be rewritten as 



)yT=^\nZ^\ = ^lnZ (M) , 



(45) 



(46) 



in terms of the variables defined in Eq. 
The free energy F of the system is 



(3F 



-lnZo^ + ^p 2 



3T 



a 2 = -\nZ ^ + I Ty\ 



(47) 



which has to be an absolute local minimum, according to Eqs. (0), in order to determine 
the equilibrium state of the system. 

The thermodynamic properties are given by 



l = -\y 2 T + lnZ ^\ H 



2„.2 



32 



T z y 



C 
k 



1 d 
~8dT 



(48) 



where S is the entropy per particle due to the orientational order, U is the energy per particle 
and C is the specific heat. They correspond to the same expressions used in Ref. |J. 

In order to study the symmetries of the partition function Z^ v ^ it is convenient to 
rewrite it as 



(Ai - A 2 ) (A x - A 3 ) (A 2 - A 3 ) 



Ai 

a: 



(9+1) A (p+i) 

-(9+1) A (P+1) 



A 



-(q+l) A (p+1) 



(49) 
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together with 

Ai = e ^ vc°s(<t> +2 r) , A 2 = e 73 W^-f ) , A 3 = e 7S (50) 
according to Eqs. ([H]), (32) and (|36|). The partition function is invariant under each 



of the following separate set of transformations, which do not change the system under 
consideration 

y^y, 0^-0, (Ai«->A 2 , A 3 -> A 3 ) , (51) 
y^y, 0^y-0, (A 3 ^A 2 , Ai->Ai), (52) 

y^y, 0^y-0, (As^Ai, A 2 ^A 2 ) (53) 

y^y, 0^0-y, (Ax -> A 3 -> A 2 -> Ax) (54) 
2tt 

y^y, 0->0 + y, (Ax^Aa^Aa^Ax). (55) 

The above set of transformations allow us to restrict the initial interval — n < < 7r to 
< < |. Furthermore, we have the following relations among the different values of A$ 

< Ax < A 2 < 1 < A 3 , < < 

o 

< Ax < 1< A 2 < A 3 , J < < (56) 

o 6 

The remaining symmetries exchange rod-like into disc-like phases 

p^q, y^y, -> vr - 0, (Ax A 2 -1 , A 3 -> A 3 _1 ) , (57) 
p^g, y^y, 0^-0, (Ai<->A 3 -1 , Aa^Ar 1 ), (58) 
p^g, y^y, 0^ + 0, (A 2 ^A 3 -\ Ax^Ax- 1 ). (59) 

The second of the above symmetries was already discussed in the classical case of Ref . || . 

For future purposes it is convenient to rewrite the ratio of the different values of Aj in 
the following way 

= e -y sin(*+f) A_2 = e _ y sin( |_0)^ Ax = sin< ^ ^ 
A 3 A 3 A 2 

7 The phases of the system 

The free energy (3F of the system, given by Eq. fl47|) , is a function of five variables; (3F = 
/3F(T,p, q, y, 0). The absolute minimum conditions (fSJ), which correspond to two equations 
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in our case, allow us to find y = y(T,p,q) and <fi = <p(T,p,q) for equilibrium, thus selecting 
the particular phase which is energetically favourable. This calculation is made numerically, 
and these results determine all the remaining thermodynamic properties of the system. The 
entropy, energy and specific heat can be subsequently calculated by using the expressions 
(|48|). The general behavior of the specific heat, as a function of temperature is presented in 
Fig. 1. Here we have labelled by T&, T u , the transition temperatures between the biaxial- 
uniaxial phase and the uniaxial-isotropic phase, respectively. Since we are only considering 
the description of the orientational modes of the system, the thermodynamic quantities 
associated with the isotropic phase will all be zero. 

7.1 The T -> (y -> oo) case. 

Let us consider the interval < <fi < | in such a way that all the ratios (|60D tend to zero 
exponentially. The corresponding limit in the partition function is taken from ( |1(J5| ), by 
observing that each of the polynomials F\ can be approximated as 

F 1 (g;A 2 ,A 3 )^A3 <? (l + ^). (61) 

In this way, the product of two F's in each term of the partition function factors out, leaving 
an additional term Fi(q; Ai, A2), which is also approximated in the way of Eq . (|B~TD . The final 
expression for the logarithm of the partition function is 

\nZj?' q) = (p + q)\nA 3 + q\nA 2 + ^ + ^ + ^. (62) 

A 3 A 3 A 2 

In order to determine which of the last three terms in the above expression dominates, we 
still have to consider the further intervals < < |, where sin0 < sin(| — 4>) < sin( : | + <fi), 
together with | < <f> < |, where sin(| — 0) < sin0 < sin(| + </>). In the first case, the free 
energy reduces to 

(3F = -(p + q)^=cos<P - g-|cos(0 - y ) - + l -Ty\ (63) 

because the term ^ dominates over the remaining fractions in (|5"2"D. Introducing the notation 



A{(j)) = 2^COS(/) + ^cos^ - ^ 



S(0) = ^sin0+^sin(0-f), (64) 

the equations that minimize the free energy are 

-A(4>) + \Ty + sine/) e"^ = 0, 

B{4>) + cos0 e- ysin(t> = 0. (65) 
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The above equations are solved by 



y = 1 ( AOo) - sin0 o e" 1 ^ 8 ^ ) . (66) 



The solution in the interval | < < | is 



</ = | f A(« - Bin(! - <M e-^-S-*)") . (67) 



This limit clearly corresponds to the biaxial equilibrium phase since A 7^ 0, (p 7^ 0. For 
the case < (j) < |, we obtain the following expression for the specific heat 

17 = 7F e T > (68) 



k VT 

where there is a mass gap given by 2q. The prediction ([68]) has to be compared with the 
classical case, where the specific heat goes like ^ = 3 + ^ s ;^)2 in that limit. Moreover, 
Eq.(|6~8"l) signals a peculiar behavior for the cases q = 0, p 7^ 0, which will be separately 
discussed in section 7.4. 

The solution ( |67T) in the interval | < < ~ leads to 

' ' 1 e"-, (69) 



k \T 

for the specific heat in the T —* limit. This is consistent with the symmetry p «-> g which 
connects the corresponding two intervals of 0. 



7.2 Biaxial-Uniaxial phase transition 

The uniaxial phase is characterized by A 7^ 0, <fi = 0. In order to determine the transition 
temperatute T b , we expand the free energy F(T, A, (p) in powers of 4> and find the correspond- 
ing minimum. The expansion is 

PF(T, A, 0) = f{z) -6z + 6Tz 2 + (z 2 g(z) + z h(z))cf) 2 + 0(</> 4 ), (70) 

where the functions f,g,h are explicitly known and z = ^. The partition function is invariant 
under the change <fi — > —</), which is reflected in the above expansion. To determine the 
transition temperature we approach from the biaxial phase and look for the corresponding 
minimum of (|70| ) having A 7^ 0, <fi 7^ 0, but small. The conditions are 

8F 

= — = 2<f)z[\g(z)+Th(z)} + 0(<J) 3 ), (71) 
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O = ^ = ^-6 + 12A + O(0 2 ). (72) 



From Eq.(frTP we obtain zg(z) + h(z) = 0, which numerical solution determines z 
Substituting this result in Eq. (|72[) the transition temperature results in 



Some results, for different values of p and q, are presented in Table I. In this table, the 
parameters \x and if) corresponding to the quantum numbers p and q are obtained from Eqs. 
(|l8l), (|19|) and (|25|) The fifth column is obtained from Eq. (ff"3|). The next columm shows 
the conversion to the classical temperature, according to Eq. ( |3"?|) . The predictions of the 
classical model are presented in the last columm. 

7.3 Uniaxial-Isotropic phase transition 

As we keep increasing the temperature we go from the uniaxial phase to the isotropic phase 
characterized by <p = 0, A = 0. To determine the transition temperature T u we proceed in 
complete analogy to the previous section. Now we expand the free energy of the uniaxial 
phase in the vicinity of z = 0, 

0F(z, 4> = 0, T) = -lnZ„ (M) + |lV =6T 2 2 - 3/ 2 z 2 - ^ I,z 3 



_ ln Q> + l)h + l)(p + g + 2) £ 108 

2 10 v 7 35 v ; 



(74) 



where we recall that z = ^. 

The extremum condition §f = requires 3Ta = dlaZ g a " ■ The general form of dlnZ ^ p ' q) 
as a function of cr is presented in Fig. 2. From here it is clear that the transition temperature 
is given by the slope of the straight line which is tangent to this curve, i.e. there is a critical 
temperature beyond which it is not possible to satisfy the consistency conditions. This 
implies the further condition 3T = ^^^r—. Both requirements imply that one is looking 
for an inflexion point of the free energy (|74"D: if = 0, = 0. 



Keeping up to fourth order terms in ([74]) , the above system is solved by 



I 1 3 I 2 

Z * = J 2 (/ 2 + 2)' T * = 2 /2 + 10 / 2 (/ 2 3 + 2)- (75) 

This solution is a good aproximation whenever z* is small. 

Nevertheless, for a range of temperatures below T*, we can verify that the uniaxial 
phase does not provide an absolute minimum for the free energy of the system. In fact, 
the calculated free energy is higher than the one corresponding to the isotropic phase, which 
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makes the uniaxial phase metastable. In this way, the true transition temperature T u between 
the uniaxial and the isotropic phase is obtained by the minimum condition together with 
the requiremente that (3F = /3F| isotropic = — ln|(p + l)(q + l){p + q + 2). In the same 
approximation as the one used in Eq. (|75|) , these new conditions provide 



4J 3 \ 4 Ij 

Zu -3I 2 (I 2 + 2Y lu ~ 2 2+ 15/ 2 (/ 2 + 2)' 

which confirm that T u < T* together with the fact that the metastable region is small indeed. 
A further discussion of the metastability region can be found in section 4.5 of Ref. || 

We can calculate the limit of the specific heat when we approach T u fom the left, in the 
case when J 3 — > (p — > q). To this end we have to find z = z(T), which we obtain by 
requiring |j = in Eq. ([74]) . We get 

/ in t \ 

(77) 




From the above equation we can verify that if we substitute the value for T u given by Eq.flTj 
we recover the value z u written in the same equation. Next we use the expression fl4"8f) for 
the specific heat to obtain 



C |t^ = 20-^- + O(/ 3 2 ). (78) 



k '— / 2 + 2 



This expression has the correct classical limit given by 9i\t^t u — 20 



The above results ([F5|), (|76|) can be further refined by taking into account the term z 5 en 
the expression ([74]). Some numerical results are presented in Table II, for different values of 
p and q. The fourth columm gives the quantum prediction for the transition temperature 
T u , which is subsequently converted to the classical temperature in the fifth columm. The 
last columm contains the classical predictions. 

7.4 The q = 0,p ^ case 

As we will see in the sequel, this situation has no counterpart in the classical description 
because here we have only an isotropic-uniaxial transition and the biaxial phase is not 
present. 

Let us consider the limit T — » within the range < <fi < | in the expression ( |103| ) for 
the exact partition function. As before, Ai, A 2 — > while A 3 —>■ oo in such a way 
that the expansion of the partition function up to second order terms leads to 



2 



-0-^.-3+ £ + 1 (£)-!(£)+.... 

The above expression differs from the q = restriction of Eq. ([62]) , in the absence of 
the previously dominating term This is because one has to take the q = value in 
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the full expression and not in a truncation of it. For example, F\(q = 0;u,v) = 1, but 
in an approximation where u > v, we can write F\ = u q (l + - + ...), which will lead to 
Fi(q — 0; u, v ) = (1 + - + . . .). It is enough to consider the expression (f79[) to first order, 
which leads to 

(3F = -J^ycoscj) - e -^ in( f-^ - e ~ ysin ^ + ^ + -Ty 2 . (80) 
V3 8 

Since we are interested in solutions near = 0, we have to keep both exponentials, which 
are of the same order now. The conditions for a minimum of the free energy are 

\Ty = -^=cos0 - sin(- - 0) e - ysin( f-^ - sin(- + 0)e~ ysin( f +*> (81) 
4 y3 3 3 

sin0 = — ( cosf- + 0)e-^5 + ^ - cos(- - (f))e- ysin ^~^) := H(<b). (82) 
p \ 3 3 / 

In writing the last equation we are assuming that ?/ ^ 0. We have that H((j) = 0) = and 
it is a direct matter to show that 4r < 0, within the interval under consideration. This 

dip ' 

means that the only solution to Eq.(|8"2"D is (ft = 0, i.e. we remain in the uniaxial phase for low 
temperatures, with no transition to the biaxial one. The remaining equilibrium condition is 

y = _-L ( p - s e -¥\ . (83) 



The above expression leads to a specific heat given by 



k=^ 6 (84) 

in the T — > 0, y — > oo limit. 

The absence of the biaxial phase in this case can also be verified from the general dis- 
cussion of section 7.2, by looking at the expansion of the free energy in powers of (p. It is 
possible to verify numerically that the coefficient of the power <fi 2 in Eq. (|70|) cannot be 
made zero, unless z = 0. That is to say, there can be no transition to the biaxial phase when 
q = 0, p^0. 

Another piece of information related to the above statement is the calculation of the 
angle if). The general result is 

p(p + 3)(2p + 3) 

cosSw = r . (85) 

2(p 2 + 3p+3)l 

Now, in the limit p —>■ oo, we recover the classical regime together with the fact that 
cos3?/> = 1, i. e. ip — 0. Recalling the phase diagram of Ref. 0, we verify that in this regime 
we have only the uniaxial and isotropic phases. 
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7.5 The p = q case 




Here I2 = p(p + 2), I 3 = 0. The specific heat, together with the entropy are shown in Figs. 
5 and 8, as a function of temperature. All the points correspond to the biaxial phase and, 
exactly as in the classical situation, the phase transition is directly to the isotropic one. In 
order to better understand the limiting values of the transition temperature, together with 
the resulting specific heats, we we will take the corresponding limit starting from the uniaxial 
phase, with J 3 = 0. In this case Eq. (fT?D, valid for any point in this phase, reduces to 

(86) 

From Eqs. (fT6|) we see that, in this limit, the uniaxial phase merges into the isotropic one, 
with z = 0, at the transition temperature T& = T u = T = y.In the conventions of Ref. [[]], 
this temperature corresponds to 

ft = 54TT" (8T) 

which reproduces the value T c — ~ for high quantum numbers. The specifc heat at the 
transition is 

k It ~ 1^+2 ~ (^TTpTT' (88) 

The above result reproduces the classical prediction. 

Let us observe that Eq. (|86|) makes sense only when T < y. In other words, we really 
approach the transition point from the sector of the uniaxial phase which has the biaxial 
phase as a boundary. 



7.6 The order parameters of the system 

In the classical case, the order parameters corresponding to the uniaxial and biaxial phases, 
respectively, are defined by [|J 

U = Tx COS ^ = 2 AcCOS( ^ b = Tx^^- = 2 A c sin( ^ ( 89 ) 
cos-0 /3J /icos-0 ' simp /3J fisinip ' 

where A c refers to the notation in this reference. The above definitions are such that U 7^ 0, 
when T <T U and 5^0 when T < T^. 
Considering that 

2\/3 2 
P= —XsiiKp a = --Acos0, (90) 

with T = jj and A c = (3JX, the corresponding order parameters in the quantum case can 
be taken as 

U=—^, 5 = --=?*—, (91) 
licosip v3/isin-0 
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in order to reproduce the classical quantities in the apropriate limit. Let us recall that fi 
and ijj are related to p and q through the equations fljjD, ([19]) and (|25|). Two complementary 
examples are presented in Figs. 10 and 11. The former corresponds practically to the 
classical situation, while the latter refers to a lower lying quantum state. 



8 Final Comments and Conclusions 

We have presented a quantum mechanical description of the orientational degrees of freedom 
of a biaxial nematic liquid in the mean field approximation, thus extending the classical 
approach proposed in Ref. ||. 

The fact that the dynamical variables carry a representation of SU(3) allows for an exact 
calculation of the quantum partition function, in terms of the Weyl formula for the characters 
of the group. In this way, the consistence equations of the mean field approximation reduce 
to two coupled algebraic equations, in a manner similar to the case in Ref. 0.The work 
presented here constitute a nice example of the applications of group theory techniques to 
quantum statistical mechanics. As expected [|13|, |14[, the classical behavior is recovered in 
the high temperature, high quantum numbers regime. In general, the discrepancies between 
the two approaches begin to appear at temperatures which are lower that the biaxial-uniaxial 
transition temperature. 

Under these circumstances and considering the actual state of the art in the production of 
liquid crystals, the liquid will probably crystalize and these effects will be irrelevant from the 
experimental point of view, in the case of the nematic phase under present-day laboratory 
conditions. This emphasizes the fact that the classical approximation is very adequate 
for the description of the phase transitions in this model. Nevertheless, further technical 
developments may lead to nematic liquids with substantially lower transition temperatures, 
such that quantum effects might turn to be important. Also, it is necessary to continue 
exploring the relevance and description of liquid crystals in high density-high temperature 
nuclear matter. 

The general phase structure of the system includes: a biaxial phase in the range < 
T < T b , a uniaxial phase in the range T b < T < T u , T b < T u , and an isotropic phase for 
T > T u . This last phase acts here as a reference, because the model does not include a 
description of translational degrees of freedom of the liquid. This structure is exhibited in 
Figs. 3, 4 and 5 (for the specfic heat) together with Figs. 7 and 9 (for the entropy). The 
specific heat always presents a discontinuity at T b , while the entropy remains continuous 
there. This corresponds to a second order phase transition. We expect the same behavior 
at the temperature T u , since both phase transitions are associated to symmetry changes in 
the liquid. The temperatures of the different phase transitions in the quantum calculation 
are about 15 — 30 % lower than the classical situation, in the most unfavourable cases (low 
quantum numbers). 

In the low temperature limit T —>■ 0, the quantum calculation of the specific heat leads to 
the result ^ ~ ^s'G T j according to Eqs. ( |68| ) and ( |69|) , signaling a mass gap characteristic 
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of the quantum spectra. The classical limit in this regime is r- = 3. Also, the quantum 
calculation is consistent with the Nernst postulate for the entropy in the T — > limit. 

Deviations from the above general picture are present in two cases: (i) the region p = q 
produces a direct transition from the biaxial to the isotropic phase, in complete analogy 
with the classical situation. This is shown in Figs. 5 and 8. (ii) the region q — 0, p ^ 
(p = 0, q 7^ 0) does not support the biaxial phase and only allows for a uniaxial-isotropic 
phase transition. This whole quantum region is the analogous to the ip = 0, (if) = |) limit 
in the classical case. These values for ij) are not allowed in the quantum regime, leading to 
q=-l(p= -1). 

Our ultimate goal, which is not discussed in this work, would be to incorporate spin 
(magnetic moment) degrees of freedom into the model for the nematic liquid, in order to 
study how the thermodynamic phases and equilibrium parameters are modified . Following 
an analogous line of development, this would require the introduction of a supergroup as 
a way of characterizing the underlying orientational degrees of freedom. In this case it 
would be also interesting to study the classical limit of the quantum version. This limit will 
shed some light into the role of the supersymmetric extension J16 of the Harish-Chandra, 
Itzykson - Zuber (HCIZ) integral, in the classical description of the model. A direct classical 
calculation of the supersymmetric model, along the lines of Ref . || , presents some difficulties, 
even though the corresponding HCIZ integral is known |Tj| . In fact, there is no Macfarlane's 
theorem in this case and the interpretation of the Grassmann numbers that naturally arise 
is not at all clear. In this situation, we expect the quantum calculation to be free of these 
ambiguities, besides of providing the complete and correct answer to the problem. 
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The Appendix 



Any set of hermitian operators S a , a = 1,2, ... ,8 satisfiying the commutators 

[S a , S b ] = if abc S c , (92) 

with f abc being the SU(3) structure constants carry a representation of the group 

The fundamental representation of SU (3) is provided by the Gell-Mann A a matrices 
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which are traceless and satisfy [A a , X b ] = if a bcK- The explicit expression for the SU(3) 
structure constants is given by 



abc 



-2itr([\ b , A c ] A Q 



(93) 



.and they are completely antisymmetric in any pair of indices. The Gell-Mann matrices also 
sat if y 

1 

{A a , A fe } = -S ab I + d abc Ac, (94) 

where d abc are completely symmetric in any pair of indices. 

The adjoint representation of SU{3) has dimension eight and possesses two invariant 
tensors 

5 ab = 2tr(X a X b ), d abc = 2tr({X b , X c }X a ). (95) 

The irreducible representations of SU(3) are labeled by two integers (p, q) where p + q 
denotes the number of boxes in the first row and q denotes the number of boxes in the second 
row, of the associated Young tableaux. The dimension D of the (p, q) representation is 



D M = ^(p + q + 2)(p+l)(q + l). 
The weights of the fundamental representation are 



1 



2'2 v / 3/ ' 



^2 



2 , 2V3J' 
20 



A*3 



o,- 



V3 r 



(96) 



(97) 



which denote de A3, Ag eigenvalues respectively. 

In the course of this research, sometimes it has proved useful to rewrite the partition 
function (fi9|) introduccing some standard polynomial functions, like Fi(k;u,v) 

k 

u k+i _ v k+i = ( u _ v ^ u j v k ~ j := {u - v)F 1 {k; u, v), (98) 

3=0 

together with F 2 (k; u,v,w), 

Fx(k] u, v) — Fi(k\ u, w) := (v — w)F 2 {k — 1; u, v, w), (99) 
which leads to the more explicit form 

k k 

F 2 (k;u,v,w) =Y f u i F 1 (k-i;v,w) = ^u^F^v, w). (100) 

i=0 i=0 

In general, for higher values of the sub index n = 3,4, . . ., the corrresponding polynomials 
can be defined inductively by 

k 

F n (k; u,v,w, . . .) = ^2^F n -i(k - i] v,w, . . .)• (101) 

8=0 

Some of the properties of the F's are the following: 

(i) They are completely symmetric with respect to the exchange of any two arguments. 

(ii) The difference between two F's which differ in just one argument is given by 

F n (k\ ill, • • • , u n , u n +i^j F n 

= («n+i - u n+2 )F n+1 (k - l;ut, . . . ,u n ,u n+1 ,u n+2 ). (102) 

In terms of the above polynomials F\ and F 2 , different cases of the partition function 
(|4"9|) are written as follows 

Zt = ° ] = F 2 (p;A u A 2 ,A 3 ). (103) 
Zt q) = F 1 (p;A 1 ,A 2 )F 1 (p;A 2 ,A3)F 1 (p;A 1 ,A 3 ) (104) 

Zt q) = x -^-(Af 1 F 1 ( ? ;A 2 ,A3)F 1 (p;A 1 ,A 3 )- 

Kf- l F x { q] A 1 ,A 3 )F l (p; A 2 , A3)) (105) 



4 P,<?) = EEA J 3 +fc (A 1 A 2 )^F 1 (p + J -^A 1 ,A 2 ). (106) 

j=0 k=0 

Some explicit expressions of the F's polynomials are the following 

Fi(l; u, v) = u + v, Fi(2; u, v) = u 2 + uv + v 2 , 

Fi(3;u,u) = u 3 + u 2 v + uv 2 + v 3 , (107) 
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F 2 (l;u,v,w) = u + v + w, 

F 2 (2; u, v, w) = u 2 + u 2 + w 2 + uv + uw + vw, 

F 2 (3;u,v,w) = u 3 + v 3 + w 3 + u 2 (v + w) +v 2 (u + w) + w 2 (u + v) 

+ uvw. (108) 
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Figure Captions 



Fig. 1 General form of the specific heat as a function of temperature, for the different 
phases. 

Fig. 2 General form of dhlZ »^ p ' q) as a function of a, for A = 0. 

Fig. 3 Specific heat for p — 17, q = 48 < — ► yu = 40.0, ■?/> = 45.0°. The upper curve is 
the classical result. The temperature is in classical units. 

Fig. 4 Specific heat for p = 25, q = 30 < — ► yu = 33.0, ip = 33.0°. The upper curve is 
the classical result. The temperature is in classical units. 

Fig. 5 Specific heat for different values of p = q < — ► ip = 30.0°, different values of ji. 
The upper curve is the classical result. The temperature is in classical units. 

Fig. 6 Specific heat for p = 1 q = 2 < — ► p = 2, 9 ip = 36.6° and p = 27, q = 41 < — ► 
/i = 40.7 ip = 36.6°. The upper curve is the classical result. The temperature is in classical 
units. 

Fig. 7 Quantum entropy for p = 25, q = 30. 

Fig. 8 Quantum entropy for p — q. The temperature is in classical units. 
Fig. 9 Quantum entropy for p = 1 q = 2 and p = 27 q = 41 . The temperature is 
in classical units. 

Fig. 10 Quantum order parameters for p = 27, q = 41. 
Fig. 11 Quantum order parameters for p — 1, q — 2. 
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Fig. 4 
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Fig. 10 
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Table I 



p 


Q 




1> 


T b 


- (T b ) c 


T b (classical) 


2 


24 


§\/709 


54.4 


2.667 


0.00282 


0.00317 


4 


43 


|V2181 


54.7 


8.000 


0.00275 


0.00287 


17 


48 


|\/3607 


45.0 


109.5 


0.02277 


0.02286 


1 


2 


fv/19 


36.6 


1.015 


0.04008 


0.05614 


27 


41 


f>/l9 


36.6 


278.3 


0.05605 


0.05614 


27 


47 


|V277 


38.6 


273.5 


0.04628 


0.04632 


25 


30 


|V2443 


32.9 


256.1 


0.07862 


0.07877 



Table II 



V 


Q 


v> 


T 


- (T u ) c 


T n (classical) 


2 


24 


54.4 


129.5 


0.1370 


0.1454 


4 


43 


54.7 


383.2 


0.1318 


0.1455 


17 


48 


45.0 


637.8 


0.1326 


0.1361 


1 


2 


36.6 


2.725 


0.1076 


0.1275 


27 


41 


36.6 


631.8 


0.1272 


0.1275 


27 


47 


38.6 


760.5 


0.1287 


0.1292 


25 


30 


32.9 


408.2 


0.1253 


0.1255 
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